Multiple-wavelength binary diffractive lenses

ABSTRACT

A dichromatic lens includes a plurality of zones being arranged on a lens structure, each of the zones having a specified radius and varying height. The lens structure focuses propagating light applicable to any intensity distribution for a plurality of wavelengths.

BACKGROUND OF THE INVENTION

The invention relates to the field diffractive optics, and in particular to multiple-wavelength binary diffractive lenses.

Diffractive optics sculpt the propagation of light to generate complex intensity and phase patterns downstream. They achieve this by imposing a particular phase and intensity pattern on the incident light. Phase-only diffractive optics, as their name implies, affect only the phase, and hence are lossless. Binary-phase diffractive optics impose only two-levels of phase. This significantly eases the fabrication of such elements. The phase shift is achieved via an optical-path difference between alternate zones. Such optics inherently exhibit chromatic aberrations.

There have been several approaches to design multiple-wavelength diffractive optics. A heterogeneous design, based on materials with differing refractive indices and dispersion to compensate for chromatic aberration, was proposed. By using phase shifts that are integer multiples of 27, harmonic diffractive lenses can be designed for specific discrete wavelengths. However, the selection of the design wavelengths is limited. A nonlinear optimization technique was used to design dual-wavelength diffractive beam-splitters. Blazed higher-order diffractive optics may also be designed for multiple wavelengths. In all these cases, the fabrication of the diffractive optic is difficult, either due to the multiple levels of phase-height or due to large aspect ratios.

SUMMARY OF THE INVENTION

According to one aspect of the invention, there is provided a dichromatic lens. The dichromatic lens includes a plurality of zones that are arranged on a lens structure, each of the zones having a specified radius and varying height. The lens structure focuses propagating light applicable to any intensity distribution for a plurality of wavelengths.

According to another aspect of the invention, there is provided a method of forming a dichromatic lens. The method includes forming a plurality of zones such that each of the zones has a specified radius. Also, the method includes varying the heights of the zones that allows focusing propagating light applicable to any intensity distribution for a plurality of wavelengths.

According to another aspect of the invention, there is provided a method of performing operations of a dichromatic lens. The method includes arranging a plurality of zones on a lens structure such that each of the zones has a specified radius and varying height. Also, the method includes lens structure focusing propagating light applicable to any intensity distribution for a plurality of wavelengths.

BRIEF DESCRIPTION OF THE DRAWINGS

FIGS. 1A-1B is schematic diagram of the optic and design methodology of the invention;

FIG. 2A is graph illustrating the transmission function of a dichromat lens structure formed in accordance with the invention, FIG. 2B is a graph illustrating the intensity in the focal plane for λ₁ illumination; FIG. 2C is a graph illustrating the intensity in the focal plane for λ₂ illumination; FIG. 2D is a graph illustrating the radial intensity distributions of the focal spots

FIGS. 3A-3B are intensity distributions in various transverse planes near the focus of a dichromat lens for λ₁ and λ₂;

FIGS. 4A-4B are graphs illustrating the standard deviations of the diffraction efficiencies at λ₁ and λ₂; FIGS. 4C-4D are graphs corresponding to data for a newly optimized dichromat lens; FIGS. 4E-4F are graphs illustrating the focal intensity distributions for the new dichromat lens at λ₁and λ₂; FIG. 4G is a graph illustrating the radial intensity distributions in the focal plane; and

FIG. 5A is a graph illustrating the transmission function of a trichromat structure formed in accordance with the invention; FIG. 5B is a graph illustrating the intensity distributions in the focal plane at the 3 wavelengths.

DETAILED DESCRIPTION OF THE INVENTION

The invention describes a technique that extends the use of nonlinear optimization to design lenses that can focus several wavelengths of light into different focal spots. The inventive lens structure focuses propagating light applicable to any intensity distribution for a plurality of wavelengths. In particular a dichromat lens is designed that focuses one wavelength, λ₁, to a bright spot and a second wavelength, λ₂, to an overlapping ring-shaped spot. The latter, with a node in its center, is a critical element in imaging schemes for breaking the far-field diffraction limit. The ring-shaped spot also has important applications in optical tweezers for trapping and manipulating particles whose refractive index is lower than that of the local environment. Such a dark spot may also have applications in trapping cold atoms. Focal spots with a dark center may be generated by focusing Laguerre-Gaussian modes, and higher-order Bessel beams.

In both cases, the null arises from the on-axis singularity in the phase of the wavefront. Such singularities can themselves be generated using diffractive elements, such as the spiral-phase plate or the spiral zone plate. The fabrication of such elements can be quite complicated, and the resulting phase profile is very sensitive to fabrication errors. Exotic interferometers have been used to generate nulls with up to 5 orders of magnitude lower intensity than the surrounding peak. These nulls are of interest in astronomy for finding faint planets orbiting a star. Phase plates that generate dark spots are also of interest in optical-projection photolithography.

One can follow the technique proposed originally by Toraldo di Francia, where the optic is composed of concentric circular zones, whose radii are the design variables. The phase shift between adjacent zones is an additional degree of freedom. This approach was shown to produce effective superresolving optical elements. Phase-only diffractive lenses with circular symmetry can be readily fabricated in a dielectric material using planar processes, enabling large arrays with high optical uniformity.

FIG. 1 is a schematic diagram of the optic and design methodology of the inventive dichromat lens structure 2. The dichromat lens structure 2 includes a multitude of zones 4 having radii r₁,r₂, . . . , r_(M) and the height of the zones, h. Stylized intensity distributions in the focal plane illustrate the design requirement for a dichromat lens that focuses a bright spot at λ₁ and a ring-shaped spot at λ₂. Outside the dichromat lens is opaque. As shown in FIG. 1, the phase shift is achieved via a varying height difference between alternate zones. The optic can be described by a circular-symmetric transmission function.

$\begin{matrix} {{T(\rho)} = \left\{ \begin{matrix} ^{\psi} & {{r_{2m} < \rho \leq r_{{2m} + 1}},{m \in I^{+}}} \\ 0 & {\rho > r_{M}} \\ 1 & {elsewhere} \end{matrix} \right.} & {{EQ}.\mspace{14mu} 1} \end{matrix}$

where ρ is the radial coordinate, r_(m), is the radius of the m^(th) zone, and M is the total number of zones. The relative phase-shift between neighboring zones, ψ, can be related to the zone height, h, via

$\begin{matrix} {{\psi = {2\pi {\frac{h}{\lambda}\left\lbrack {{{Re}\left( {n(\lambda)} \right)} - 1} \right\rbrack}}},} & {{EQ}.\mspace{14mu} 2} \end{matrix}$

where Re(n(λ)) is the real part of the refractive index of the lens material. The dichromat lens can focus propagating light applicable to any intensity distribution for a plurality of wavelengths beside 2 described above.

The Fresnel-Kirchoff formulation of the scalar-diffraction problem to model the propagation of light from the optic to the plane of observation is used. A normally incident uniform plane normally incident uniform plane wave is assumed for simplicity. The intensity in the observation plane is then given by

$\begin{matrix} {{{\mathrm{\Upsilon}\left( {\rho^{\prime},z,\lambda,\left\{ {\overset{\rightarrow}{r},h} \right\}} \right)} = \left| {\frac{1}{i\; \lambda}{\int_{0}^{\infty}{\rho \ {\rho}{\int_{0}^{2\pi}\ {{\varphi}\; {T\left( {\rho,\left\{ {\overset{\rightarrow}{r},h} \right\}} \right)}\frac{\exp \left( {\frac{2\pi}{\lambda}\sqrt{\rho^{2} + \rho^{\prime 2} - {2{\rho\rho}^{\prime}{\cos (\varphi)}} + z^{2}}} \right)}{\sqrt{\rho^{2} + \rho^{\prime 2} - {2{\rho\rho}^{\prime}{\cos (\varphi)}} + z^{2}}}\left( {1 + \frac{z}{\sqrt{\rho^{2} + \rho^{\prime 2} - {2{\rho\rho}^{\prime}{\cos (\varphi)}} + z^{2}}}} \right)}}}}} \right|^{2}},} & {{EQ}.\mspace{14mu} 3} \end{matrix}$

where z is the propagation distance along the optical axis, and ρ′ and φ′ are cylindrical coordinates in the observation plane. The design variables are the radii of the zones, {right arrow over (r)}={r₁,r₂, . . .,r_(m)} and the height of the zones, h.

The key step of the design process is the nonlinear optimization. The goal of the optimization is to achieve a certain diffraction pattern in the focal, or observation, plane. This goal is described in terms of an energy function. The technique is illustrated via a dichromat lens that focuses λ₁ to a round spot and λ₂ to a ring-shaped spot. The energy function for this design is then expressed as:

E({{right arrow over (r)},h})=−ω₁∫_(o) ^(ρ′) ¹ Υ(ρ′,λ₁ ,{{right arrow over (r)},h})ρ′dρ′−ω ₂∫_(ρ′) ₂ ^(p′) ⁸ Υ(ρ′,λ₂ ,{{right arrow over (r)},h})ρ′dρ′+ω ₃ Υ(ρ′=0,λ₂ ,{{right arrow over (r)},h}),   EQ. 4

where ρ₁′ is the nominal radius of the round spot at λ₁, and ρ₂′ and ρ₃′ are the nominal inner and outer radii of the ring-shaped spot at λ₂. The variables w1, w2 and w3 are positive weights that allow relative emphasis of one term or another in the energy function. The last term adjusts the depth of the null in the center of the λ₂ spot. The optimization algorithm attempts to minimize the energy function. In addition, the following constraints are imposed.

τ_(p)>τ_(q)∀M≧p>q>1, {p,q}∈ I   EQ. 5

τ_(p)−τ_(p−1)>Δ>0∀M≧p>1, p ∈ 1   EQ. 6

τ₁>2Δ  EQ. 7

The first constraint ensures that the zones are retained in the correct order during optimization. The second constraint ensures that the width of each zone is greater than Δ, a constraint dictated by fabrication technology.

This inventive technique can be extended to describe the case of oblique illumination, as shown in FIG. 1B. In many optical applications, it is important to design the lenses with large acceptance angles. This is achieved by simply using a different set of equations to compute the intensity distribution in the observation or focal plane. All other aspects of the technique remain the same.

For an oblique angle of incidence, α, the intensity in a plane at a distance z from the lens is given by

$\begin{matrix} {{{Y_{\alpha}\left( {\rho^{\prime},z,\lambda,\alpha,\left\{ {\overset{\rightarrow}{r},h} \right\}} \right)} = {\frac{A}{i\; \lambda}{\int_{0}^{2\pi}\ {{\varphi}{\int_{0}^{\infty}{\rho^{\prime}\ {\rho^{\prime}}{T\left( {\rho,\left\{ {\overset{\rightarrow}{r},h} \right\}} \right)}\frac{\exp \left( {\; k} \right)}{}\left\{ {{\cos \; \alpha} + \frac{z}{d}} \right\}}}}}}},} & {{EQ}.\mspace{14mu} 8} \end{matrix}$

where T is the transmittance of the lens, (ρ,φ) are radial co-ordinates in the observation plane, and (ρ,φ) are radial co-ordinates in the lens plane, and d is given by

$\begin{matrix} {d = \sqrt{\left( {{\rho \; \cos \; \varphi} - {\rho^{\prime}\cos \; \varphi^{\prime}}} \right)^{2} + \left( {{\rho \; \sin \; \varphi} - {\rho^{\prime}\sin \; \varphi^{\prime}}} \right)^{2} + z^{2}}} & {{EQ}.\mspace{14mu} 9} \end{matrix}$

If an incoherent addition of angles is desired, the energy function can be modified to sum over a range of incident angles. In a general form, the energy function can then be written as:

$\begin{matrix} {{E = {- {\sum\limits_{i}^{\;}{w_{i}{\sum\limits_{l}^{\;}{\int_{r_{i\; 1}}^{r_{i\; 2}}{{Y_{\alpha}\left( {\rho^{\prime},z,\lambda,\alpha,\left\{ {\overset{\rightarrow}{r},h} \right\}} \right)}\rho^{\prime}{\rho^{\prime}}}}}}}}},} & {{EQ}.\mspace{11mu} 10} \end{matrix}$

where the negative sign indicates that the energy is being minimized, w_(i) are the weights associated with wavelength, λ_(i), r_(1i) and r₂₁ are the inner and outer radii of the ring into which λ_(i) is focused (if that is what is desired), and the summation over ₁ takes the range of incident angles into account.

Although, the scalar Fresnel-Kirchoff diffraction theory is used in the examples, the invention is equally applicable for any theory that models the propagation of light from the diffractive lens to the observation or focal plane. Other applicable theories are the first and second Sommerfeld diffraction equations, finite-difference-time-domain methods, or the like. non-radial lenses as well. Also, the same design technique can be applied to linear (one-dimensional) lenses for 1-D focusing. The inventive technique can be applicable to any intensity distribution for each wavelength that can be specified by the user.

The transmission function of the optic is described by a piece-wise linear function as seen in equation (1). Free-space propagation is a highly nonlinear function of the spatial co-ordinates. Furthermore, the energy function and the constraints add additional nonlinearity. These characteristics make this problem ideal for a genetic algorithm. The genetic algorithm is an iterative mathematical version of natural selection. At each iteration, individuals from the population are chosen to mate based on the values of their energy functions. Offspring are produced by sharing “genes”, i.e., the variables {{right arrow over (r)},h} or mutation, i.e., a random perturbation of the variables {{right arrow over (r)},h}, This procedure is repeated until a set of variables is found that gives a global minimum for the energy function. The genetic algorithm is particularly appropriate to solve problems that are not well suited for standard optimization algorithms, such as when the energy function is discontinuous, non-differentiable, or highly nonlinear.

The technique described above was used to design a dichromat lens with numerical aperture 0.7, focal length 100 λ₁, and 80 zones. The design wavelengths were λ₁=400 nm, and λ₂=532 nm. The optimization was carried out with ρ₁′=ρ₂′=0.5 λ₁/NA and ρ₃′=1.22 λ₂/NA-ρ₂′. The lens was assumed to be made of polymethylmethacrylate (PMMA), however other similar materials can be used. The refractive indices of PMMA were measured as 1.501 (λ₁=400 nm) and 1.487 (λ₂=532 nm). Δ was set to 0.5 λ₁. The optimum weights were empirically determined, and set as w₁=1, w₂=10 and w₃=10⁴.

The transmission function of the dichromat lens after 25 generations of a genetic algorithm is shown in FIG. 2A. Within the aperture of the dichromat lens, the dark rings are phase shifted with respect to the white ones; their relative height difference being 0.7367 λ₁. Outside the dichromat lens is opaque. The intensity distributions in the focal plane are shown in FIGS. 2B and 2C for λ₁ and λ₂, respectively. Radial cross-sections through the centers of the focal spots are plotted in FIG. 2D. Clearly, the technique was successful in designing a dichromat lens that focuses λ₁ to a round, bright spot and λ₂ to a ring-shaped spot. The inset in FIG. 2D shows that the intensity in the central null is about 5 orders of magnitude below that in the surrounding peak.

The focusing characteristics of the dichromat lens change as the observation plane is swept through the focus. This defines the useful depth-of-focus of the lens. FIG. 3A shows the focused spots for various values 30-38 of z. The ring at λ₂ is well-defined for a range of approximately λ₁/2. FIG. 3B shows the corresponding data 40-48 for λ₁. It was noticed that the intensity in the central node at λ₂ increases rapidly with defocus.

Assuming that the dichromat lens would be fabricated using planar processes, it is important to understand the sensitivity of the focusing characteristics to errors introduced during fabrication. Fabrication errors manifest themselves as errors in the radii of the zones and the height of the zones. Their effect can be simulated by adding randomly generated errors to the zone-radii and the phase height.

{right arrow over (r)}={r ₁ ,r ₂ , . . . ,τ _(M) }→{right arrow over (r)}={r ₁+{tilde over (δ)}_(r1) ,r ₂+{tilde over (δ)}_(r2) , . . . ,r _(M)+{tilde over (δ)}_(rM)}  EQ. 11

h→ h=h+ δ _(h)   EQ. 12

where {tilde over (δ)}_(rp) and {tilde over (δ)}_(h) are randomly generated from two normal distributions of zero mean, and standard deviations, δ_(r) and δ_(h) respectively. In order to quantify the effect of the error, one can calculate the distribution of the focusing efficiencies of the error-prone dichromat lenses at the two wavelengths, and then, calculated their corresponding standard deviations. The focusing efficiencies are defined as:

$\begin{matrix} {{{\eta_{1}\left( \left\{ {\overset{\rightarrow}{r},h} \right\} \right)} = \frac{\int_{0}^{\rho_{1}^{\prime}}{{\mathrm{\Upsilon}\left( {\rho^{\prime},\lambda_{1},\left\{ {\overset{\rightarrow}{r},h} \right\}} \right)}\rho^{\prime}{\rho^{\prime}}}}{\int_{0}^{\infty}{{\mathrm{\Upsilon}\left( {\rho^{\prime},\lambda_{1},\left\{ {\overset{\rightarrow}{r},h} \right\}} \right)}\rho^{\prime}{\rho^{\prime}}}}}{and}} & {{EQ}.\mspace{14mu} 13} \\ {{\eta_{2}\left( \left\{ {\overset{\rightarrow}{r},h} \right\} \right)} = \frac{\int_{\rho_{2}^{\prime}}^{\rho_{3}^{\prime}}{{\mathrm{\Upsilon}\left( {\rho^{\prime},\lambda_{2},\left\{ {\overset{\rightarrow}{r},h} \right\}} \right)}\rho^{\prime}{\rho^{\prime}}}}{\int_{0}^{\infty}{{\mathrm{\Upsilon}\left( {\rho^{\prime},\lambda_{2},\left\{ {\overset{\rightarrow}{r},h} \right\}} \right)}\rho^{\prime}{\rho^{\prime}}}}} & {{EQ}.\mspace{14mu} 14} \end{matrix}$

and The robustness of the dichromat lens was investigated by calculating the standard deviations of the diffraction efficiencies, σ{η₁({{right arrow over ({tilde over (r)},{tilde over (h)}})} and σ{η₂({{right arrow over ({tilde over (r)},{tilde over (h)}})} as a function of δ_(r) and δ_(h). The results for the original dichromat lens are shown in FIGS. 4A and 4B. In order to design a more robust dichromat lens, one can define a new energy function as follows:

E _(r)({{right arrow over (r)},h})=μ{E({{right arrow over ({tilde over (r)},{tilde over (h)}})}+σ{E({{right arrow over ({tilde over (r)},{tilde over (h)}})},   EQ. 15

where μ{ } and σ{ } denote mean and standard deviation respectively. This energy function was used to design a new dichromat lens with the same parameters as described earlier. Standard deviations, δ_(r)=λ₁/8 and δ_(h)=λ₁/16 were used during the optimization. After 25 generations of the genetic algorithm, a second dichromat lens design was obtained. This design was significantly more robust to fabrication errors as illustrated by the dramatically reduced variation in their diffractive efficiencies, as shown in FIGS. 4C and 4D. The focal intensities of the new dichromat lens at λ₁ and λ₂ are shown in FIGS. 4E and 4F. Note that the depth of the central null at λ₂ is slightly worse than in the original dichromat lens. This is characteristic of optimization, where one quality is traded off against another. It is likely that better energy functions can attain higher quality nulls, while maintaining error-tolerance. FIG. 4G is a graph illustrating the radial intensity distributions in the focal plane. The inset 60 shows the λ₂ distribution on a log-scale.

The proposed technique is easily extended to more than two wavelengths. One can demonstrate this by designing a trichromat, a lens that can focus λ₁ and λ₃ into bright spots, while λ₂ is focused into a ring-shaped spot. The energy function is modified as follows.

E({{right arrow over (r)},h})=−w ₁∫_(o) ^(ρ′) ¹ Υ(ρ′,λ₁,{{right arrow over (r)},h})ρ′dρ′−w ₂·_(p′) ₂ ^(ρ′) ³ Υ(ρ′,λ₂ , {{right arrow over (r)},h})ρ′dρ′+w ₂ Υ(ρ′=0, λ₂ ,{{right arrow over (r)},h})−w ₄ ∫_(o) ^(ρ′) ⁴ Υ(ρ′,λ₃ , {{right arrow over (r)},h})ρ′dρ′,   EQ. 16

where ρ₄′ is the desired spot radius for λ₃ and w₄ is the weight for the last term. The parameters for the optimization were λ₃=633 nm, ρ₄′=1.22 λ₃/NA and w₄=−10. The refractive index of PMMA at λ₃ was measured as 1.4812. All other parameters were the same as for the original dichromat lens.

FIG. 5A shows the transmission function of the trichromat after 50 generations of optimization. FIG. 5B shows the intensities in the focal plane for the three wavelengths. The inventive technique can, therefore be extended to an arbitrary number of discrete wavelengths with the added flexibility of being able to specify different focal distributions at each design wavelength.

The invention describes a technique based on non-linear optimization using genetic algorithms to design binary, phase-only diffractive optics for multiple wavelengths. It has been demonstrated the efficacy of this technique by designing a dichromat lens, a lens that focuses λ₁ to a central bright spot and λ₂ to a ring-shaped spot with a deep central null. It is shown that the design technique is flexible enough to incorporate robustness to fabrication errors, and is easily extended to more than two wavelengths simply by incorporating appropriate optimization criteria.

Although the present invention has been shown and described with respect to several preferred embodiments thereof, various changes, omissions and additions to the form and detail thereof, may be made therein, without departing from the spirit and scope of the invention. 

1. A dichromatic lens comprising: a plurality of zones arranged on a lens structure, each of said zones having a specified radius and varying height, said lens structure focusing propagating light applicable to any intensity distribution for a plurality of wavelengths.
 2. The dichromatic lens of claim 1, wherein said lens structure focusing propagating light to a bright spot at a specified first wavelength and a ring-shaped spot at a specified second wavelength when a phase shift occurs because of varying heights of said zones.
 3. The dichromatic lens of claim 1, wherein the radius and height of each of the zones is based on the desired intensity distributions in the observation plane for any number of said wavelengths.
 4. The dichromatic lens of claim 1, wherein lens structure receives incoming light at an oblique angle.
 5. The dichromatic lens of claim 1, wherein said lens structure comprises polymethylmethacrylate (PMMA).
 6. The dichromatic lens of claim 1, wherein said lens structure comprises an aperture.
 7. The dichromatic lens of claim 2, wherein said ring-shaped spot comprises dark rings that are phase shifted with respect to white rings using said aperture.
 8. The dichromatic lens of claim 2, wherein said ring-shaped spot comprises a deep central null.
 9. The dichromatic lens of claim 2, wherein said first wavelength is UV and said second wavelength is the visible wavelength.
 10. A method of forming a dichromatic lens comprising: providing a lens structure; forming a plurality of zones on said lens structure such that each of said zones has a specified radius; varying the heights of said zones that allows focusing propagating light applicable to any intensity distribution for a plurality of wavelengths.
 11. The method of claim 10, wherein said lens structure focusing propagating light to a bright spot at a specified first wavelength and a ring-shaped spot at a specified second wavelength when a phase shift occurs because of varying heights of said zones.
 12. The method of claim 10, wherein the radius and height of each of the zones is based on the desired intensity distributions in the observation plane for any number of said wavelengths.
 13. The method of claim 10, wherein lens structure receives incoming light at an oblique angle.
 14. The method of claim 10, wherein said lens structure comprises polymethylmethacrylate (PMMA).
 15. The method of claim 10, wherein said lens structure comprises an aperture.
 16. The method of claim 11, wherein said ring-shaped spot comprises dark rings that are phase shifted with respect to white rings using said aperture.
 17. The method of claim 11, wherein said ring-shaped spot comprises a deep central null.
 18. The method of claim 11, wherein said first wavelength is UV and said second wavelength is the visible wavelength.
 19. A method of performing operations of a dichromatic lens comprising: arranging a plurality of zones on a lens structure such that each of said zones having a specified radius and varying height; focusing propagating light applicable to any intensity distribution for a plurality of wavelengths.
 20. The method of claim 19, wherein said lens structure focusing propagating light to a bright spot at a specified first wavelength and a ring-shaped spot at a specified second wavelength when a phase shift occurs because of varying heights of said zones.
 21. The method of claim 10, wherein the radius and height of each of the zones is based on the desired intensity distributions in the observation plane for any number of said wavelengths.
 22. The method of claim 10, wherein lens structure receives incoming light at an oblique angle.
 23. The method of claim 10, wherein said lens structure comprises polymethylmethacrylate (PMMA).
 24. The method of claim 10, wherein said lens structure comprises an aperture.
 25. The method of claim 11, wherein said ring-shaped spot comprises dark rings that are phase shifted with respect to white rings using said aperture.
 26. The method of claim 11, wherein said ring-shaped spot comprises a deep central null.
 27. The method of claim 11, wherein said first wavelength is UV and said second wavelength is the visible wavelength. 